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We study the semileptonic B s and B decays into the 2317) and Dq( 2400) resonances, respec¬ 
tively. With the help of a chiral unitarity model in coupled channels we compute the ratio of the 
decay widths of both processes. Using current values of the width for the B° —> Dq(2400) + i>i D we 
make predictions for the rate of the B° —> D* 0 (2317) + vi l~ decay and for the DK invariant mass 
distribution in the B° —» D K vil~ decay. 

PACS numbers: 


I. INTRODUCTION 

The recent discovery of many mesons with charm contributed to the revival of hadron spectroscopy (for recent 
reviews, see Ref. [1]). Two interesting examples of these mesons are the D* 0 (2317) [2] and Dg(2400) [3] scalar 
resonances. As it happened to other states, their measured masses and widths do not match the predictions from 
potential-based quark models. This disagreement motivated several non-conventional (exotic) interpretation of these 
states. Among them the most popular are multiquark configurations in the form of tetraquarks or meson molecules 
[1]. Since their masses are located below the DK and D S K thresholds it is quite natural to think that they are bound 
states of DK and D S K meson pairs. In the case of the D* 0 ( 2317), additional support to the molecular interpretation 
came recently [4] from lattice QCD simulations. In all previous lattice studies of the D* 0 (2317), it was treated as a 
conventional quark-antiquark state and no states with the correct mass (below the DK threshold) were found. In Ref. 
[4], with the introduction of DK meson operators, the right mass was obtained. In [5] the scattering length of KD 
from QCD lattice simulations was extrapolated to physical pion masses and then, using the Weinberg compositeness 
condition [6, 7] the dominance of the KD component in the D* 0 (2317) state was concluded. A reanalysis of the results 
of [4] has been done in [8] using the information of the three energy levels of [4] and going beyond the effective range 
formula. The dominance of the DK component of the D* 0 (2317) was firmly established in [8]. On the other hand the 
analysis of the D* 0 (2317) -7 D *7 radiative decay suggests that the D* 0 (2317) is most likely an ordinary cs state [9], 
although this is disputed in [10]. It is therefore important to carry out further studies to clarify this question. One 
aspect to be investigated is the production of D* 0 (2317) and Dg(2400) in B s and B decays, respectively. 

The dominant decay channel of the B s meson is into the D s meson plus anything. Therefore various important 
properties of the cs mesons can be studied in the B s weak decays. In particular, they can shed more light on the 
controversial D* 0 (2317) meson, whose nature is still under debate, as discussed above. In recent years there has been a 
significant experimental progress in the study of the properties of the B s mesons. The Belle Collaboration considerably 
increased the number of observed B s mesons and their decays [11]. Moreover, B s mesons are copiously produced at 
Large Hadron Collider (LHC) and precise data on their properties have been taken by the LHCb Collaboration [12]. 
New data are expected in near future [13]. The study of weak B s decays is primarily devoted to the improvement 
in the determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, but there are several interesting 
topics in hadron physics to be investigated in these processes. 

The D* 0 (2317) has been measured mostly in B factories and probably because of its narrow width (T < 3.8 MeV) 
it has not been seen in some channels. One of them is the B® semileptonic decay, i.e., B° s —> D*q“ v 1 1~ (Fig. 1). There 
are several theoretical estimates of the branching fraction of this decay channel [14-19] and they predict numbers 
which differ by up to a factor two. In all these calculations, a source of uncertainty is in the hadronization. In the 
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present work, assuming that the £>2 0 (2317) and £>2(2400) are dynamically generated resonances, we try to improve 
the hadronization and the treatment of final state meson-meson interactions. Moreover, in order to further reduce 
uncertainties we compute the ratio of decay widths: 


R = 


S»^D; 0 +(2317) I 7 i i- 




ui /- 


(1) 


We calculate the left side of this equation and then, using the available experimental information about the process 
B° £>2(2400) + Hi l~ (see Fig. 2), we can extract ^ b °^ d * (2317)+ ■ We consider also the B~ —> £12(2400)° vi l~ 

decay (Fig. 3). The formalism is very similar to the one presented in Ref. [20, 21] for nonleptonic B decays. 

As it is depicted in Fig. 4, after the W emission the remaining c — q pair is allowed to hadronize into a pair of 
pseudoscalar mesons (the relative weights of the different pairs of mesons is known). Once the meson pairs are produced 
they are allowed to interact in the way described by chiral unitarity model in coupled channels and automatically the 
£>2 0 (2317) and £>2(2400) resonances are produced. 



B° 

S 


—--«i--— 

s s 

FIG. 1: Semileptonic decay of into Vil~ and a primary cs pair. 



*4 


d d 

FIG. 2: Semileptonic decay of B° into vil~ and a primary cd pair. 

This paper is organized as follows. In Sec. II, we formulate the semileptonic B decay widths into D resonances and 
give our model of the hadronization. Next in Sec. Ill we consider D resonance production via meson coalescence after 
rescattering, and in Sec. IV we calculate the production of two pseudoscalars with prompt production plus rescattering 
through a D resonance. Then in Sec. V we formulate meson-meson scattering amplitudes to generate the £>2o(2317) 
and £>2(2400) resonances. In Sec. VI we show our numerical results of the semileptonic B decay widths. Section VII 
is devoted to drawing the conclusion of this study. 
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FIG. 3: Semileptonic decay of B into vil and a primary cu pair. 


II. SEMILEPTONIC B DECAYS 

Let us first formulate the semileptonic B decays into D resonances in the following decay modes: 

B° s -+ £>* 0 (2317)+^r, 

B° —>■ Dq(2A00) + DiI~ , (2) 

B~ Do(2400)°F ; r, 

where the lepton flavor l can be e and g. For this purpose we express the decay amplitudes and widths in a general 
form in Sec. IIA and give our model of the hadronization in Sec. IIB. 


A. Semileptonic decay widths 


In general, by using the propagation of the W boson and its couplings to leptons and quarks, we can express the 
decay amplitude of B —> Dl~ hadron(s), Tg, in the following manner: 


■rr - ■ 9w a 1 - 75 
- iT b =uh-j=j —-— v v x 


^9cx.f3 


P 2 - 


— .gwVbc a 1 — 75 / - Tx 

X Ucl ^ 2 ^ ' ^^had)? 


( 3 ) 


where it;, v v , it c , and Ub are Dirac spinors corresponding to the lepton l - , neutrino, charm quark, and bottom quark, 
respectively, gw is the coupling constant of the weak interaction, Vb c is the Cabibbo-Kobayashi-Maskawa matrix 
element, and Mw is the W boson mass. The factor T4ad consists of the wave function of quarks inside the B meson 
and the hadronization contribution in the final state, and it will be evaluated in the sections below. In the following 
we neglect the squared momentum of the W boson (p 2 ) which should be much smaller than M^ in the B decay 
process, and therefore the decay amplitude becomes 


■ j' _ Vbc x cu r\ 

J- B — ^ y/y ^ o 


X Vhadj 


( 4 ) 


where we have introduced the Fermi coupling constant G F = and defined the lepton and quark parts 

of the W boson couplings as. 


L a = Ul"f a (l - 7 5 )lV, Qol = Ucla (1 - 75 )ub, 


( 5 ) 


respectively. 

Now let us calculate the decay widths of the semileptonic B mesons into D resonances. In the calculation of the 
decay widths, we take the absolute value of the decay amplitude Tg and average (sum) the polarizations of the 
initial-state quarks (final-state leptons and quarks). Therefore, in terms of the amplitude in Eq. (4), we can obtain 
the squared decay amplitude as 


5 S' 1 '* i 2 


|G F I4 c 14ad| 2 




pol 


4 


( 6 ) 
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where the factor 1/2 comes from the average of the bottom quark polarization. Then, by using the conventions of 
the Dirac spinors and traces of Dirac matrices summarized in Appendix A, we can calculate the lepton part of the 
amplitude (5), which reads 


Y L a L t/3 =tr 

pol 


an n , \ pfa + mi 

7 (1 - 75) ^r (1+7s)7 


PuPi + P?Pt - Pv ■ Pi9 ap - ie 


=2 


a/3 _ n f P a °P 


PvpPlo 


m^rrii 


(7) 


where p u and pi (m v and mi) are momenta (masses) of the neutrino and lepton l , respectively. In a similar manner, 
we can calculate the quark part of the amplitude, which is given by Eq. (5) 


Y Q*Ql =tr 

pol 


7a (1 -75) 


j) b + m b j) c + m c 


2 m b 


-(1 + 75)7/3- 


2 m c 


P nrCT 


=2 


PbaPcft “ 1 “ Pca.Pb(3 Pb ' PcQa.(3 ^^potcrfiPbP. 
m b m c 


( 8 ) 


with the momenta (masses) of the bottom and charm quarks, p b and p c ( m b and m c ), respectively. Now we take 
a heavy quark limit and assume that the momentum of both the bottom and charm quarks are zero at the B rest 
frame. Then we have p b = ( m b , 0 ) and p% = (m c , 0 ), and the quark part of the W coupling can be rewritten as 


^ ' QaQp — 25a[ 3 ? 

pol 


(9) 


at the B rest frame. Here we note that this is the delta function rather than the metric g a p. As a consequence, the 
square of L a Q a with polarization summation gives 


Y^ \L a Q a \ 2 =4 

pol 


(28 a pp“pf - p„ ■ PlS a pg a P) B rest 


m„mi 


16 (E v Ei)b rest 


m v mi 


where we have used 5 a pg a P = g 00 + g 11 + g 22 + g 33 = —2. Finally we obtain the squared decay amplitude: 


(I>bI 2 = 


pol 


4|G ? Fl4cl'had|" 

m v m\ 


(E v E t ) 


B rest• 


( 10 ) 


(ii) 


With the above squared amplitude we can compute the decay width. We will be interested in two types of decays: 
three-body decays, such as B 3 —> Df 0 vi l ~, and four-body decays, such as B 3 —> D + K° Ui l~ and also for the similar 
B° and B~ initiated processes. As it will be seen, both decay types can be described by the amplitude Tb with 
different assumptions for I4ad- The final formulas for 3 and 4-body decays are then given by: 

r 3 = 


r 4 = 


respectively. In the equations, ms, m„, and mi are respectively the masses of the B meson, neutrino and lepton l, 
T 3 ( 4 ) is the three- (four-) body decay amplitude, and the summation symbols represent the average of the polarizations 
in the initial state and the sum over the polarizations in the final state. Moreover, the n-body phase space d$> n has 
been introduced as 


4m v mi 
Zms J 

[ d * 3 ££ i^ai 2 , 

(12) 

pol pol 


^m„mi 

2m B J 

f d^YYM 2 ’ 

pol pol 

(13) 


n 

- n 


d 3 pi 1 

(2tt)3 2 El 


(2tt) 4 6 4 (pb — ptot), 


(14) 
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where Ei = ■\/p‘f + mf is the on-shell energy of i-th particle with its mass p^ is the four-momentum of the initial 
B meson, and p^ ot is the sum of the final-state momentum: 

n 

Pto t = Yl p i ’ Pi = t Ei ' ( 15 ) 

2=1 

In order to proceed with the calculation we need a prescription of hadronization, i.e., after the W emission in Figs. 
1-3 we must specify a way to convert the outgoing quarks into hadrons and compute I4ad- This will be done in the 
next subsection. 


B. Hadronization 

The conversion of quarks into hadrons in the final stage of hadron reactions is a long-standing problem which 
up to now has no definitive solution. Since the energies involved are of the order of a few GeV or less, this is a 
non-perturbative process. For particles produced in very high energy collisions and with high transverse momentum, 
we can use fragmentation functions, which are extracted from data phenomenologically and then refined with a 
perturbative QCD treatment. In some case one can develop an approach based on effective Lagrangians [22], In the 
process considered here, in contrast to the high energy case where many particles are produced along with the formed 
hadron, only one quark-antiquark pair is produced and hadronization is mostly a recombination process which binds 
together the existing quarks. Here we follow [20] and describe hadronization as depicted in Fig. 4. An extra qq 
pair with the quantum numbers of the vacuum, uu + dd + ss + cc, is added to the already existing quark pair. The 
probability of producing the pair is assumed to be given by a number which is the same for all light flavors and which 
will cancel out when taking ratios of decay widths. We can write this cq (uu + dd + ss + cc) combination in terms of 



cq (uu+d d + ss+cc) 


FIG. 4: Schematic representation of the hadronization cq —> cq (uu + dd + ss + cc). 


pairs of mesons. For this purpose we follow the work of [21] and define the qq matrix M: 

( uu ud us uc \ 

du dd ds dc 

su sd ss sc 

y cu cd cs cc J 


M = 


which has the property 

M ■ M = M x (uu + dd + ss + cc). 
Now, in terms of mesons the matrix M corresponds to [23] 

/>° 


7^+7^' 

7T + 

K+ 

D° \ 

7r 

_L^o + J_ ?7 + _L„' 

K° 

D~ 

K~ 

K° 

CS 

+ 

/ip' d 7 

D° 

D+ 

Df 

Pc / 


V 

Hence, in terms of two pseudoscalars we have the correspondence: 


V3 


(16) 


(17) 


(18) 


cs (uu + dd + ss + cc) = (<f> ■ cj )) 43 = D°K + + D + K° + Df 


jp' j +VcDf 


(19) 
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FIG. 5: Diagrammatic representation of D*q (2317) production via meson coalescence after rescattering. 


+ -W) + D+K° + Vc D+ (20) 

cu ( uu + dd + ss + cc) = (</> ■ />) 41 = D° ( — -=n° 4 — -=r] -\ — -=rf ) + D + n~ + D+K~ + rj c D° (21) 

V02 03 06 ) 

for £>* 0 (2317) + , Uq( 2400) + , and Uq( 2400)° production, respectively. Then, for simplicity we concentrate on the 
relevant channels for the description of the D resonances. In fact, it was pointed out in Ref. [24] that the most 
important channels for the description of D* 0 (2317) (Dq( 2400)) are DK and D s rj (Dn and D S K). Therefore, the 
weights of the channels to generate the D resonances can be written in terms of the ket vectors as 

|(#) 4 3 > = V2\DK(0 , 0)) - ^\DsV(0, 0)), (22) 

1(00)42} = -^fliMl/2, 1/2)} + \D s K(l/2, 1/2)), (23) 

I(00)4i) = yf|lMl/2, -1/2)) - |iWl/2, -1/2)), (24) 

where we have used two-body states in the isospin basis, which are specified as (/, I3) and are summarized in 
Appendix A. We note that, due to the isospin symmetry, both the charged and neutral Dq (2400) are produced with 
the weight of | (</></>) 42 ) = — | (</></>) 4 i), which means that the ratio of the decay widths into the charged and neutral 
Dp (2400) is almost unity. By using these weights, we can express Vhad in terms of two pseudoscalars. 

Once the quark-antiquark pair hadronizes into two mesons they start to interact and the D resonances can be 
formed as a result of complex two-body interactions with coupled channels described by the Bethe-Salpeter equation. 
If the resonance is formed, independent of how it decays, the process is usually called “coalescence” [25] and it is 
a reaction with three particles in the final state (see Fig. 5). If we look for a specific two meson final channel we 
can have it by “prompt” or direct production (first diagram of Fig. 6 ), and by rescattering, generating the resonance 
(second diagram of Fig. 6 ). This process is usually called “rescattering” and it is a reaction with four particles in the 
final state. Coalescence and rescattering will be discussed in the next sections. 


cd (uu + dd + ss + cc) = ( (j> ■ 4>) i2 = D°tt + + D + (-p 7 r° + 


02 03 


III. COALESCENCE 

In this section we consider D resonance production via meson coalescence after rescattering (see Fig. 5). This 
process has a three-body final state with a lepton, its neutrino and the resonance R. The hadronization factor, Vhad: 
can be obtained as 


Wad(A>: 0 (2317)) = C ^ 0 2Gdk9dk — -j^GD sV gD a 


(25) 







FIG. 6: Diagrammatic representation of DK production: directly (on the left) and via rescattering (on the right) in B° decays. 


I4ad (£>5(2400)+) 


-4iad(-D5(2400)°) = C 


~~\l -^GnngD-K + G DsK g DsK 


(26) 


Here y, is the coupling constant of the D resonance to the i -th two meson channel and Gi is the loop function of two 
meson propagators (see Sec. V) 


Gi(s) = i 


d 4 q 1 1 

(2tt) 4 q 2 — to 2 + it (P — q) 2 — to' 2 + it 


(27) 


where P^ is the total four-momentum of the two-meson system, and thus P 2 = s with s being the invariant mass 
squared of the two-meson system, and rrii and to' are the masses of the two mesons in channel i. An important point 
to note is that the prefactor C is the same in all decay modes and contains dynamical factors common to all reactions 
only because we are assuming that in the hadronization the SU(3) flavor symmetry is reasonable, i.e., the quark pairs 
cs (Fig. 1) and cd (Fig. 2) hadronize in the same way. We further assume that C is a constant and therefore is 
canceled when we take the ratio of decay widths as in Eq. (1). 

Now we can evaluate the decay widths by using the formula of Eq. (12). Inserting Eq. (25) [or Eq. (26)] into Eq. (11) 
and the latter into Eq. (12) we can write the decay width of the D resonance production via meson coalescence (Fig. 5) 
as 


lf;oal — 


TO„TO; 

1287t 5 to^ 




dfln / dH, 


4|G F HcWad(J*)P 

m u mi 


(£,£)) 


B rest 


(28) 


where p^ is the momentum of the D resonance in the B rest frame and p u is the momentum of the neutrino in the 
vl rest frame, both of which are evaluated as 


« cm — 
Pd ~ 


a 1/2 K b . wS'lt O 

2 m B 


(29) 


Pv = 


ml 


2 M 


(vl) 



(30) 


with the Kallen function A(x, y, z) = x 2 + y 2 + z 2 — 2 xy — 2 yz — 2 zx and the D resonance mass m#. The tilde on 
characters indicates that they are evaluated in the vl rest frame unless explicitly mentioned. The solid angles Qd and 
are for the D resonance in the B rest frame and for the neutrino in the vl rest frame, respectively, and is 

the vl invariant mass. After performing the angular integrals, we obtain the final expression of the decay widths for 
the coalescence of the D resonance: 


Tcoal — 


| G F l 4 c l 4 a d (£>*)| 2 r &l) - 

0 o o / ^^inv Pd Pv\P J i'P J l)B rest? 


2n 3 m 2 b 


(31) 


where the integral range of is [mi + m v , ms — mu]. In the equation, (E„Ei)b rest is the product of the v and l 

energies averaged over the neutrino solid angle and it is calculated in the following way. Before the angular integral, 
we have an exact relation: 


(E l/ Ei'j b res t 


( Pv-Pb)(pi ■PB ) 


(32) 

















Since ( p v ■ pb) and (pi ■ pb) are Lorentz invariant, we may evaluate them in the vl rest frame as 


Pv,i ■ Pb = E v jE b - pu,i - Pb, Pv = ~Ph 
For simplicity we neglect lepton masses, so we have 


E v = Ei = p v = pi = 


Ml 


("0 


ipu,i = \Pv,i\)- 


(33) 


(34) 


On the other hand, the kinetic condition leads to an exact form of Eb- 


+ [M^] 2 - m 2 R 

PjB — — 


2 Ml 


(vl) 


In this way we have 


Then (E„Ei)b rest becomes 


Pv,i ■ Pb = 


1% + lML 


*B ~r l M inv } 2 - m R 


Pv Pb 


Pv,i ■Pb ■ 


‘b + [ML 


Pv • Pb 


(EvE t ) 


B rest — 


m B + i M inv } 2 “ (P" ' PB? 


\ 4tob ) m| j 

Integrating the second term with the neutrino scattering angle in the vl rest frame, we have 


1 


d cos 9, 


(. PvPb ) 2 1 (PvPbY 


i-i 


3 m 2 B 


where ps = y E B — m 2 B . As a result, we obtain ( E„Ei)b rest as 


b + [ A 4S 2 -n4 \ 1 ( PvPb) 2 


(EvE t ) 


B rest — 


4 tob 


3 m 2 B 


(35) 


(36) 


(37) 


(38) 


(39) 


IV. RESCATTERING 


Next, the production of two pseudoscalars with prompt production plus rescattering through a D resonance is 
calculated with the diagrams shown in Fig. 6, and its hadronization amplitude Vhad in the isospin basis is given by 


bhad(-DA') = C ^V2 + V2GdkTdk^dk — Gd sV Td srj^DK^j , 
14ad (D s rj) = C + VzGdkTdk^Dsti — -^jMGd^Tds^Ds^ , 


(40) 


(41) 


V had (Z?7T) =C\- 


o ~~ \J 2 + G DaK T DaK ^ Dn j , 


(42) 


Vh a d (D s K) — C [ 1 — \j -G D7T T Dn ^ DaK + G DsK T DsK ^ DaK 


(43) 
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Again we see that the prefactor C is the same in all the reactions. In order to calculate decay widths in the particle 
basis, we have to multiply by the appropriate Clebsch-Gordan coefficients. 

Inserting Eq. (40) [or Eqs. (41), (42), and (43)] into Eq. (11) and the latter into Eq. (13) we can derive the differential 
decay width dTi/dM^, where i represents the two pseudoscalar states and M^ v is the invariant mass of the two 
pseudoscalars, as 


dTi |G F I4 c Vhad(i)| s 


dM 


(<) 


87r 5 m^ 


dM^P cm p u pi(E u E t ) B rest , 


(44) 


where P cm is the momentum of the vl system in the B rest frame, p v is defined in Eq. (30), and pi is the relative 
momentum of the two pseudoscalars in their rest frame, both of which are evaluated as 


pcm _ 


("0l 
inv 

2 m B 


[m£>] 2 ) 


(45) 


Pi = 


A 1 / 2 ([M ; 1 


inv 1 2 > m l 


2 M. 


(0 



Here we note that pi is a quantity in the rest frame of the two pseudoscalars rather than the vl system. 


(46) 


V. THE DK-Dst] AND Dir-D S K SCATTERING AMPLITUDES 

In this section we will discuss in more detail the amplitudes which appear in Eqs. (40), (41), (42), and (43). We 
formulate meson-meson scattering amplitudes for the rescatterings to generate the _D* 0 (2317) and Hq( 2400) resonances 
in the final state of the B decay. In Ref. [24] it was found that the couplings to DK and D s rj are dominant for U* 0 (2317) 
and the couplings to Dir and D S K are dominant for Dq( 2400). Therefore, in the following we concentrate on DK-D s r\ 
two-channel scattering in isospin 1 = 0 and Dir-D S K two-channel scattering in I = 1/2, extracting essential portions 
from Ref. [24] and assuming isospin symmetry. Namely, we obtain these amplitudes by solving a coupled-channel 
scattering equation in an algebraic form 


Tij(8) = Vi 3 (s) + Y, V ik (s)G k {s)T kj (s), (47) 

k 

where i, j. and k are channel indices, s is the Mandelstam variable of the scattering, V is the interaction kernel, and 
G is the two-body loop function. 

The interaction kernel V corresponds to the tree-level transition amplitudes obtained from phenomenological La- 
grangians developed in Ref. [24]. Here we summarize the tree-level amplitude in the isospin basis (for the two-body 
states in the isospin basis, see Appendix A). Namely, for the DK-D s rj scattering in / = 0 we have 

V DK n DK( s ’ u ) = - q 7-7 b(t -u)+s-u + m 2 D + m\\ , (48) 

'JjnJD 


KKYni 8 , u ) = KdTkd( S , t,u) = - 


phen. 


6V3 Uf D 


[7 (u — t) - (3 + 7 )(s - u) - m 2 D - 3 m 2 K + 2m 2 ], 


(49) 


V dIT'dsv^ s ' = ~ Qf-fo + 2t - u ) + 2m D + 6m K - 4ml\ , 


(50) 


and for the Dir-D S K scattering in I = 1/2 we have 

1 


V Dl e Dn( s ’ u ) = - l2 f J D + (7 + 4)(s — u) + 2m 2 D + 2m 2 ], 


U phen ' -(s t u) = UP nen - ( s t u) = 


phen. 


2 VoUf D 


[ 7(1 — u) + s — u + m 2 D + m 2 K ] . 


(51) 
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TABLE I: Pole position ^/Spoie, coupling constant (/*, compositeness Xi, and elementariness Z for the D resonances in the 
isospin basis. 


Dio 

(2317) 

Hq (2400) 

^pole 

2317 MeV 

^pole 

2128 - 160i MeV 

9dk 

10.58 GeV 

9Dn 

9.00-6.18i GeV 

9Dsv 

-6.11 GeV 

9d s k 

-7.68 + 4.35z GeV 

Xdk 

0.69 

X £>71- 

0.34 + 0.41* 

Xd s ti 

0.09 

X d s k 

0.03 - 0.12* 

Z 

0.22 

Z 

0.63 - 0.28* 


t t phen• 

V D S K D S K 


(s, t, u ) = 


1 

6 firfD 


[ 7 (t — u) + s- u + m 2 D + 2 m 2 K - ml] . 


(52) 


where t and u are Mandelstam variables. In these equations, /„■ and fp represent the pion and D meson decay 
constants, respectively, and m w , mx , and mu are the masses of pion, kaon, and D mesons, respectively. In addition, 
in order to treat effectively interactions of heavy mesons, we have introduced a parameter 7 = (wl/wh ) 2 as the 
squared ratio of the masses of the light to heavy vector mesons (respectively tol and mu), which are exchanged 
between two pseudoscalar mesons. Then we perform the on-shell factorization and the s-wave projection to give the 
interaction kernel V in Eq. (47): 


1 7 1 

V(s) = - / dcos9 V phen '(s, t(s, cos9), u(s, cos9)), 

2 7-i 


(53) 


where 9 is the scattering angle in the center-of-mass frame. 

For the loop function G, on the other hand, we use the expression in Eq. (27). In this study we employ the 
dimensional regularization, so we can express the loop function as 


Gfc(s) = 


167T 2 


Ofc(A<: 


reg 


. ml 

) + —~ 

/i r eg 


s + m fc - m k 
2 s 


In ■ 


"k 


2A 1 / 2 ( i 




Si artanh (Vi 


l k> 


mi 


+ m 12 — s 


, (54) 


with the regularization scale /i reg and the subtraction constant a*,, which becomes a model parameter. In this approach, 
D resonances can appear as poles of the scattering amplitude T,j(s) with the residue gigj'. 


Tij(s) 


9i9j 

S Spole 


+ (regular at s = s po i e )- 


(55) 


The pole is described by its position s po i e and the constant gi, which can be interpreted as the coupling constant 
of the D resonance to the i channel. In this study only the subtraction constant in each channel is the model 
parameter. Actually, the meson masses are fixed as m n = 138.04 MeV, mx = 495.67 MeV, m v = 547.85 MeV, 
mu = 1867.23 MeV, and mo, = 1968.30 MeV, and we take 


U = 93 MeV, f D = 165 MeV, 


(56) 


m L = 800 MeV, m H = 2050 MeV, (57) 

for the pion and D decay constants and masses of the light and heavy vector mesons, respectively. On the other 
hand, the subtraction constant, as a model parameter, is determined so as to generate a pole of H* 0 (2317) at the 
right place, i.e., to reproduce the mass reported by the Particle Data Group from the square root of the pole position, 
-Spoie- In this study we assume that all the subtraction constants take the same value for simplicity, and employ 
anx = od.ij = cld-k = a D s K = —1-27 at /x reg = 1500 MeV. Indeed, with these values of the subtraction constant we 
obtain the pole positions listed in Table I. The values of the coupling constants g. t are also given in Table I. 

As one can see from Table I, the D* 0 (2317) state has zero decay width with Im^/s po i e = 0, since we do not include 
the D s 7 T° decay channel, for which the isospin symmetry breaking is necessary. As a result, the coupling constants 
also become real and have positive values. The DK coupling constant is about two times larger than that of the D s rj 
coupling constant, and their values are in agreement with the results obtained in Ref. [24]. On the other hand, the 
Dq( 2400) state has a decay width (2 Im^/Spoie ~ 320 MeV) to the Dir decay channel, and the coupling constant is 
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TABLE II: Ratios of decay widths and branching fractions of semileptonic B decays. 


R 

0.45 

Lfl- -tD* (2400)° Vi l - /Lb°->£>* (2400)+ Pg- 

1.00 

B[B° -»• £>3(2400)%/“] 

3.0 x 10“ 3 (input) 

B[B~ ->• £>3(2400)%/“] 

3.2 x 10“ 3 


complex. The magnitude of the Dtt coupling constant is larger than that of the D S K coupling constant, and they are 
very close to the values in Ref. [24]. In the following we will use the coupling constants of the D resonances in Table I 
for the coalescence of D resonances in the semileptonic B decays [see Eq. (31)] and use the scattering amplitude for 
the meson-meson invariant mass distributions of the D resonances [see Eq. (44)]. 

Let us further discuss the structure of the D resonances in this model from the point of view of compositeness, 
which is defined as the contribution from the two-body part to the normalization of the total wave function and 
measures the fraction of the two-body state [26-30]. Actually, the coupling constant gi is found to be the coefficient 
of the two-body wave function in Refs. [31, 32], and the expression of the compositeness in the present model is 


Xi = -at 


dGj 

ds 


(58) 


J ® — ^pole 

On the other hand, the elementariness Z, which measures the fraction of missing channels, are expressed as 


Z 


E 


9j9i 


G^G, 
ds 


-J S — Spole 


(59) 


We note that in general both the compositeness Xi and elementariness Z become complex values for a resonance 
state and hence one cannot interpret the compositeness (elementariness) as the probability to observe a two-body 
(missing-channel) component inside the resonance. However, a striking property is that the sum of them coincides 
with the normalization of the total wave function for the resonance and is exactly unity: 


Y,Xi + Z= 1, (60) 


which is guaranteed by a generalized Ward identity proved in Ref. [33]. Therefore one can deduce the structure by 
comparing the value of the compositeness with unity, on the basis of the similarity to the stable bound state case. 
The values of the compositeness and elementariness of the D resonances in this approach are also listed in Table I. 
The result indicates that the £>* 0 (2317) resonance, which is obtained as a bound state in the present model, is indeed 
dominated by the DK component. This has been corroborated in the recent analysis of QCD lattice results of [8]. In 
contrast, we may interpret that the £>3(2400) resonance is constructed with missing channels, although the imaginary 
part for each component is not negligible. 


VI. NUMERICAL RESULTS 

Now we show our numerical results of the semileptonic B decay widths. As we have seen, we fix the hadronization 
process of the two mesons in Sec. IIB and we employ an effective model in Sec. V so as to determine the strength of 
the couplings of the D resonances to the meson-meson channels. In this way, we can calculate the ratio of the decay 
widths in the coalescence treatment as well as in the rescattering. 

First we consider the coalescence case. The numerical results are summarized in Table II. The most interesting 
quantity is the ratio R = Tg 0 ^ D * ( 23 i 7 )+p,;- /^'b°^d*( 2400 )+DiI~ i n the coalescence treatment, which removes the 
unknown factor C in the hadronization process. The decay width in the coalescence is expressed in Eq. (31). The 
coupling constants of the two mesons to the D resonances are determined in Sec. V and listed in Table I. We emphasize 
that we have no fitting parameters for the ratio R in this scheme. As a result, we obtain the ratio of the decay widths 
as R = 0.45. On the other hand, we find that the ratio Tb-^d*( 2 A 00 )°viI~ /r b°^>d*( 2400 )+& t i- is 1-00, which can be 
expected from the same strength of the decay amplitude to the charged and neutral £>3(2400) due to the isospin 
symmetry, as discussed after Eq. (24). 

Then, we can fix the absolute value of the common prefactor C by using experimental data of the decay width. 
Actually, the branching fraction of the semileptonic decay B° —> £>3(2400)"%U to the total decay is reported as 
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TABLE III: Branching fraction of the process B° — > DJ 0 (2317) + uil in percentage 


Approach 

B[B° S -+ D: 0 (2317)+h i r] 

This work 

0.13 

QCDSR + HQET [14] 

0.09 - 0.20 

QCDSR (SVZ) [15] 

0.10 

LCSR [17] 

0.23 ±0.11 

CQM [16] 

0.49 - 0.57 

CQM [18] 

0.44 

CQM [19] 

0.39 


(3.0 ± 1.2) x 10 -3 by the Particle Data Group [34]. By using this mean value we find C = 7.28, and the fractions 
of decays B® — > D* 0 (2317) + iil~ and B~ — > Dq (2400) o P;Z _ to the total decay widths are obtained as 1.3 x 10 -3 and 
3.2 x 1CV 3 , respectively. The values of these fractions are similar to each other. The difference of the fractions of 
B° — > Dq( 2400) + vil~ and B~ —> Dq(2400) 0 7'zU comes from the fact that the total decay widths of B° and B~ are 
different. 

In Table III we compare our predictions for B[B° S —>■ D* s0 {2317) + vil~] with the results obtained with other ap¬ 
proaches. Although not explicitly mentioned by the authors, from the reading of the works we can see that we should 
attach an uncertainty of at least 10 % to the numbers without theoretical error bars. The discrepancy between the 
calculated branching fractions can be of a factor five, showing that there is a large room for improvement on the 
theoretical side. Our approach is the only one where the _D* 0 (2317) + is treated as a mesonic molecule. Looking at 
Table III we can divide the results in two groups: the first four numbers, which are “small” and the last three, which 
are “large”. In the second group, the constituent quark models (CQM) yield larger branching fractions. Understand¬ 
ing the origin of the discrepancies requires a very careful comparative analysis of all the ingredients of the different 
approaches and it is beyond the scope of the present work. However it is tempting, as a first speculation, to attribute 
these differences to the differences in spatial configurations, which are inherent to each approach. The parent B s 
meson is a compact state, with a typical radius of the order of the lowest charmonium radius, i.e. (r) — 0.4 fm. In 
constituent quark models all the D s mesons, being relatively heavy cq states, should also be compact and hence the 
overlap between the initial and final state spatial wave functions is large. In the molecular picture of the _D* 0 (2317) + , 
a bound state of two mesons is expected to have a large radius, of the order of a few fm, and therefore in this case 
the overlap between initial and final of wave functions is small, reducing the corresponding branching fraction. In the 
QCD sum rules formalism, there is no explicit mention to the spatial configuration of the interpolating currents and 
it is difficult to say anything. 

Next we consider the rescattering process for the final-state two mesons formulated in Sec. IV. We use the common 
prefactor C = 7.28 fixed from the experimental value of the width of the semilcptonic decay B° —> Dg(2400) + 7'//”. 
The meson-meson scattering amplitude is obtained in Sec. V, and we further introduce the D s n° channel as the 
isospin-breaking decay mode of D* 0 (2317). Namely, we calculate the scattering amplitude involving the D s 7 r° channel 
as 


T, 


i—>D s 7T° 


9i9D s n° 


s-[M D , o -iT D .j2Y 


(61) 


for i = DK and D s rj. We take the Z3* 0 (2317) mass as Md* = 2317 MeV, while we assume its decay width as 
T d*q = 3.8 MeV, which is the upper limit from experiments [34]. The D* 0 (2317)-« coupling constant gi (i = DK, 
D s r]) is taken from Table I, and the D* 0 (2317)-D s 7r° coupling constant gn^ 0 is calculated from the Z?* 0 (2317) decay 
width as 


9d 3 7T° 



(62) 


with the pion center-of-mass momentum p w , and we obtain gu s Tv ° = 1-32 GeV. 

The results of the differential decay width dTi/dM^ (44), where i represents the two pseudoscalar states, are 
shown in Fig. 7. The figure is plotted in the isospin basis. Therefore, when translating into the particle basis we use 
the relation according to the weight of states given in Appendix A: 

[D°K+] = [D+K°} = l -[DI <], 


( 63 ) 
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2 2.2 2.4 2.6 2.8 3 

M^l [GeV] 

FIG. 7: Differential decay width dTi/dM^l for the two pseudoscalars channel i in the isospin basis. Here we consider the 
semileptonic decays B° —> (DK) + uil ~, (D a n°) + vil~ and B° —> {Dn) + vil~ . The DK and D s n° channels couple to the 

-Ds 0 (2317) + resonance, and Dn to the Dq( 2400) resonance. The peak height for the D s n° channel is dr 0s7r o /dM^ an ~ 10 -13 . 


[G+tt 0 ] = [D s n°], 

(64) 

[D°n + } = 2[D + n°} = ^[Dn], 

(65) 


where [ AB] is the partial decay width to the AB channel. An interesting point is that the DK mode shows a rapid 
increase from its threshold « 2360 MeV due to the existence of the bound state, i.e., the D* 0 (2317) resonance. In 
experiments, such a rapid increase from the DK threshold would support the interpretation of the D* 0 (2317) resonance 
as a DK bound state. The strength of the DK contribution in the M^ v > 2.4 GeV region is similar to that of Dir, 
which corresponds to the “tail” for the Dq( 2400) resonance. In fact, the position of the Gq( 2400) peak in Fig. 7 might 
be shifted to higher invariant masses, since we have underestimated the mass of the Dg(2400) resonance in our model 
compared to the experimental values 2318 MeV and 2403 MeV for neutral and charged Gq( 2400), respectively [34] 
(note that the experimental uncertainties in the position and width of this resonance are large). On the other hand, 
the D s ir° peak coming from the G* 0 (2317) resonance is very sharp due to its narrow width. The D s ir° peak height 
is about 30 times larger than the Dn one coming from Hq( 2400), but when integrating the bump structure of the 
differential decay widths we obtain a ratio of semileptonic B decays into D* 0 (2317) to Gq( 2400) close to 0.45, as 
obtained in the coalescence treatment above. 

The spectra shown in Fig. 7 are our predictions and they may be measured at the LHCb. They were obtained in 
the framework of the chiral unitarity approach in coupled channels and their experimental observation would give 
support to the G* 0 (2317) and Dq( 2400) as dynamically generated resonances, which is inherent to this approach. 


VII. CONCLUSION 

We have extended the formalism developed in [20] and applied it to semileptonic B and B s decays into resonances, 
which are interpreted as dynamically generated resonances. As in [20], we start studying the weak process at the quark 
level and, as a “final state interaction”, the outgoing quark-antiquark pair couples to meson pairs, which rescatter and 
form resonances, which then decay in well defined channels. This process, discussed in Sec. IIB, is a very economic 
hadronization mechanism with a single parameter, C. After fixing it with the help of experimental information on 
the B° —> Uq( 2400) + P;U decay, we make predictions for the semileptonic decay width of the G* 0 (2317), shown in 
Table III and also for the invariant mass spectra shown in Fig. 7. 

We have added a new information related to the nature of the G* 0 (2317) as an object with a dominant DK molecular 
component, which is the DK mass distribution in the B® —> ( DK) + vil~ decay. The simultaneous measurement of 
the decay rate into the D* 0 (2317) resonance and the related DK mass distribution are hence strongly encouraged 
to gain further knowledge on the nature of this resonance. The experimental confirmation of our predictions would 
give additional support to the D* 0 (2317) and .Dp(2400) resonances as dynamically generated resonances from the 
meson-meson interaction. 
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Appendix A: Conventions 


In this Appendix we show conventions used in this study. Throughout this article we employ the metric in four¬ 
dimensional Minkowski space defined as g = gp v = diag(l, —1, —1, —1) and the Einstein summation convention is 
used unless explicitly mentioned. 

We introduce the Dirac spinors u(p , s) and v(p, s), where p is three-momentum of the field and s represents its 
spin, as the positive and negative energy solutions of the Dirac equation, respectively: 


(j/> — m)u(p , s) =0, {j/> + m)v(p, s) = 0. 


(Al) 


Here m is the mass of the field, j) = YPn with being the Dirac gamma matrices, and p p ^ \Jp 1 + m 2 , p'j is the 

on-shell four-momentum of the solution. In this study the Dirac spinors are normalized as follows: 


u(p, s)u(p, s') = 5 ss i, v(p, s)v(p, s') = -S ss >, 
with u = idy 0 and v = idy 0 , and hence we have 


Y u (P’ s ) m (P’ s ) = 2m . 

s 

Y v ( p ’ s ) = ^ 7 "- 


The trace identities used in this study is summarized as follows: 

tr [WW] - 4 (g^g»° - + g^g^), 


(A2) 


(A3) 


(A4) 


tr [7 5 7 M 7 v 7 P 7 a ] = -4 ie^ pa , 


(A5) 


t r = t r YlsYYY] = 0, 

where 75 = ij 0 "/ 1 "/ 2 ! 3 an( 4 e M1/pcr is the Levi-Civita symbol with the normalization e 0123 = 1. 

The phase convention for mesons in terms of the isospin states | J, I 3 ) used in this study is given by 

k + ) = -|l, 1), \K~) = 11/2, —1/2), 

I D°) = 11/2, -1/2), 


(A6) 


(A7) 


while other meson states in this study are represented without phase factors. As a result, we can translate the 
two-body states used in this study into the isospin basis, which we specify as (/, I3), as 


\DK(0, 0)) = -7=1 D°K+) + Y\ D +K°), 


(A8) 


\DsV(0, 0)) = \Df v ), 


(A9) 
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\Dn(l/2, 1/2)} = -^| \D°n+) + -^=|^ + 7r°), (AlO) 

\D s K(l/2, 1/2)) = | DfK°), (All) 

\Dn(l/2, —1/2)) = +^\D°n°) + /fp+O, (A12) 

\D s K(1/2, -1/2)) = -| D+K~). (A13) 

At last we summarize the Feynman rules for the weak interaction used in this study. We express the Wvl coupling 
as 


- jgw ^ 1-75 

lv Wvl “ 4 rxi o ) 


V 2 


(A14) 


with gw being the coupling constant of the weak interaction and the Wbc coupling as 




■ gwVbc a 


1-75 

2 


(A15) 


where Vb c is the Cabibbo-Kobayashi-Maskawa matrix elements, whose absolute value is |14 c | ~ 0.041. The W boson 
propagator with four-momentum is written as 


= p^-M^+,0 ' <A16) 

with the mass of the W boson Mw- The coupling constant gw and the mass of the W boson Mw are related to the 
Fermi coupling constant G f as 


G f 


gw 

4V2M2, 


1.166 x 1CT 5 GeV“ 2 . 


(A17) 
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